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Abstract 

Inflation driven by a single, minimally coupled, slowly rolling field generically yields a 
negligible primordial non-Gaussianity. We discuss two distinct mechanisms by which a non- 
trivial potential can generate large non-Gaussianities. Firstly, if the infiaton traverses a 
feature in the potential, or if the inflationary phase is short enough so that initial transient 
contributions to the background dynamics have not been erased, modes near horizon-crossing 
can acquire significant non-Gaussianities. Secondly, potentials with small-scale structure 
may induce significant non-Gaussianities while the relevant modes are deep inside the hori- 
zon. The first case includes the "step" potential we previously analyzed while the second 
"resonance" case is novel. We derive analytic approximations for the 3-point terms gener- 
ated by both mechanisms written as products of functions of the three individual momenta, 
permitting the use of efficient analysis algorithms. Finally, we present a significantly im- 
proved approach to regularizing and numerically evaluating the integrals that contribute to 
the 3-point function. 
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1 Introduction 

The primordial power spectrum or 2-point correlation function of the temperature anisot- 
ropies in the Cosmic Microwave Background (CMB) is well-measured out to large multipoles. 
The higher moments of a distribution are, in general, independent of the 2-point function, 
but the CMB anisotropics are at least approximately Gaussian. Theoretically, we know that 
non-Gaussian component to the CMB will always be induced by non-linear gravitational 
couplings between modes after they reenter the horizon pj while single field slow-roll infiation 
yields a primordial non-Gaussian component roughly 100 times smaller than that induced 
by gravitational mode-couplings [^i3J. The two terms are additive, so recovering this latter 
primordial signal is next to impossible, even before contending with cosmic variance. Multi- 
field models may generate larger non-Gaussianities [H El [6l [71 [8] , but this is by no means 
a generic property of these scenarios. Consequently, the detection of a primordial 3-point 
function would immediately falsify a very large class of infiationary models. Conversely, non- 
slow-roll models with higher order derivative terms, such as DBI infiation P, [101 [HI [12] and 
k-infiation [131 [13!, typically generate large non-Gaussianities [15]. Further references, 
include more complicated multi-field, non-local or ghost theories, can be found in Refs. [TH 
[IIl[l8l[T9l[20l[211[22l[23l[2l]. In this paper we investigate simple models - in the sense that 
they are driven by a single, minimally coupled scalar field with a canonical kinetic term - 
which generate substantial non-Gaussianities. 

Constraining the non-Gaussian signal in a CMB dataset is a highly non-trivial problem. 
Firstly, it depends on the choice of estimators. At the moment, only two concrete estimators 
have been constructed: the fj^f and /X' forms [251 |26l [23 [281 [291 [301 [31] , and both are 
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scale- invariant. The recent claim of a detection of a non-zero 3-point function [32] in the 
WMAP 3-year data ^28j relies on the estimator developed in Ref. [211 ESI ES] which is of 
the "local" form [331 [25]. However non-Gaussianities that have strong scale dependence are 
not well-described by /j^™' and /^^"\ and will require scale-dependent estimators, based 
on theoretically motivated ansatzen for the primordial non-Gaussianities. In addition to 
computing the primordial 3-point term, comparing the predictions of a specific inflationary 
model to the CMB requires us to evolve this signal through to recombination, and to project 
it onto the sphere on the sky, which is a convolution of the primordial 3-point term with 
3 /-valued spherical Bessel functions. In general, this process is computationally expensive 
but simplifies dramatically if the 3-point function has special algebraic properties [3^ [35]. 

Using Maldacena's elegant formalism [3], the primordial 3-point curvature correlation 
function is computed via a set of integrals (over time) of products of three mode functions 
(or their derivatives) and slow roll parameters. Looking at these integrals, we identify two 
mechanisms which can create substantial non-Gaussianities. The first class consists of po- 
tentials with a localized violation of slow roll. This can take the form of a small localized 
feature, including the step models [Ml [37] whose 3-point term was first accurately computed 
by the present authors in [38]. In addition models with a short inflationary phase can have 
initial transients in their dynamics, which will not be fully erased before the longest modes 
leave the horizon. In these cases, the 3-point term for modes which are leaving the horizon 
during the violation of slow roll can be magnified by three orders of magnitude, without 
ruining the fit to the 2-point function. The second class arises when a small ripple is super- 
imposed on top of an otherwise smooth potential. This induces a "resonance" inside one of 
these integrals, giving the 3-point function a substantial amplitude before the modes cross 
the horizon. In the former case the 2-point function may be substantially modified; in the 
second case the modification of the 2-point function is very small, even though the 3-point 
function is large. The latter mechanism has not been described previously and the seemingly 
contrived field theoretic potential may in fact arise naturally in brane inflation models [35] . 

With this information in hand, we construct rough analytic approximations to the corre- 
sponding 3-point terms. These expressions have the factorizable form required by [M], which 
means that we can efficiently compute constraints on step potentials or similar models from 
CMB data, although at this point we are only interested in a qualitative match to the numer- 
ically computed 3-point term. We check our semi-analytic estimates for the 3-point function 
using an improved version of the code described in [38j, which has a much cleaner scheme 
for removing numerical divergences in the 3-point integrals. The integrands are products of 
large, rapidly oscillating terms. Analytically, these are finite, but their numerical evalua- 
tion is non-trivial, and we show how to transform them into an explicitly finite form before 
the numerical evaluation is carried out. From a practical perspective, this means we can 
compute the 3-point function for "triangles" which contain two very different scales. 

The paper is organized as follows. In Section [2] we review the computation of the 3-point 
function, and in Section [3] we describe the numerical algorithm. Sections [Hand [5] discuss the 
two distinct mechanisms for generating large primordial non-Gaussianities within minimally 
coupled single scalar field infiation, and we derive the approximate analytic ansatzen for 
the 3-point function. We then use our numerical methods to compute the exact 3-point 
function for these models, and show that these approximations yield fair representations of 
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the underlying non-Gaussianities. In Section [H] we summarize and discuss our results and 
plans for future work. We follow the notational conventions of ^38j and set the reduced 
Planck mass Mp to unity except when presenting final results or defining parameter values. 



2 Single scalar field inflation models and Maldacena's formalism 

Consider scalar field inflation with an arbitrary potential 



S 



During inflation, spacetime is described by the Friedman-Robertson- Walker metric 
ds^ = -dt^ + a^{t){dr'^ + r'^dn) = a\-dT^ + dr^ + r^dQ) , 



(2.1) 



(2.2) 



and the conformal time r runs from — oo to 0. Dots denote derivatives by cosmic time t 
while primes denote derivatives with respect to r. 

The evolution of the single scalar field is described by 



a 



1 dV 



0" + 2-0' + -— = 

a dcp 

while the scale factor obeys the Friedmann equation 

'1 



(2.3) 
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(2.4) 



The solution of the scalar field to these coupled set of equations is the trajectory of the 
scalar field in phase space (0(t),0(t)). We can also describe this trajectory with the slow- 
roll parameters [40j 
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(2.5) 
(2.6) 



Slow-roll infiaton occurs when |e| -C 1 and \ri\ <^ 1. 

Since we are interested in models with complicated potentials, we will need to solve 
their perturbation equations numerically, before deriving analytical approximations to the 
3-point functions. Moreover, we will need to track the perturbations while they are inside 
the horizon, and thus cannot resort to the large wavelength approximation [8[ El HJl H3] . 
Instead we use the approach introduced by Maldacena [3] , in what follows we use the notation 
of both this paper and Ref . [55] . 

We first split the Hamiltonian into its quadratic component Hq and its higher order 
interaction component i//, 



(2.7) 
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The superscript / signifies that these modes are evolved using the hnear (i.e. free field) 
equations of motion. Since we are interested only in the 3-point correlation functions, Hi 
contains terms up to third order in linear variables. 

One can then show (see the Appendix of [H]) that any given combination of product of 
fields Q, can be evolved by a simple unitary transform 



Qit) 



Texp 



I I Hj{t)dt 

to 



Q\t) 



Texp -i / Hi{t)dt 



to 



(2. 



Here T and T refer to anti-time-ordering and time-ordering, but since we have only third 
order terms, the time ordering is not a factor in the equation. Using Eq. (12.81) [3J, the 3-point 
correlation function for the Bardeen curvature C [35] is 



(2.9) 



(C(r, ki)C(r, k2)C(r, kg)) = -i / dr' a ([C(r, ki)C(r, k2)C(r, kg), Hi{t')]) 



TO 



This is well-suited to numerical calculations, as we now simply need to solve for the Fourier 
mode of the linear order perturbation in order to compute the relevant integrals. 

The main advantage of Maldacena's approach is the usage of the ADM formalism, where 
the constraint equations can be conveniently solved and interaction terms expanded. The 
computation of Hi for a minimally coupled single field inflationary model was performed in 
detail in Ref. [3]. It has the general form 



(2.10) 



where ^ denotes either C^' or 5^, and gi are functions of the slow-roll parameters e{t) 
and Ti{t). We see that the coupling strength of the interaction Hamiltonian depends on the 
dynamics of the background encoded by the slow-roll parameters. We emphasize that this 
equation is exact to all orders in slow-roll parameters, and does not depend on the slow-roll 
conditions. In other words, even if we violate the slow roll conditions these equations still 
hold. In order to compute the 3-point function, we must evaluate the set of integrals listed 
in Appendix \^ It is now clear why standard single field slow roll inflation does not generate 
large non-Gaussianities: the interaction couplings are functions of the slow-roll parameters, 
and hence are small by construction. However, large non-Gaussianities are possible if these 
couplings behave in a non-trivial way while keeping the viability of the power spectrum. 



3 Numerical integration 

Our goal is to integrate Eqs. (12.91) and (I2.10p numerically. We first numerically solve the 
background and linear order equations of motion for the fields using the free field Hamiltonian 
- to do this we solve the relevant equations of motion using a standard Runge-Kutta 6th order 
integrator. We then plug the solution into Eq. (12.91) and integrate them from — oo < r < 
to obtain the 3-point function. 
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The integrals (see Appendix |X] for the exphcit formulas) possess the generic form 



/ oc 3? 



'^e.nd 



Wui{Tend) j dTa'^g{e,r],r]')^i{T)^2iT)^3{T) 



(3.1i; 



where ^ can denote ({k) or d({k) / dr, and g is some function of the slow roll parameters which 
differs from term to term. These integrals are formally convergent in the limit r — > — cxd, 
but cutting them off at a finite value of r exposes the oscillatory nature of the integrand 
whose amplitude blows up rapidly as r grows large and negative, introducing a spurious 
contribution of 0{1). Physically, when the modes are well within the horizon, they oscillate 
rapidly compared to the rate of change of the interaction terms, so the contribution to the 
integrals almost cancels. We can see this by rotating the integrals r r(l + ie) into the 
imaginary plane, giving the oscillatory part of the integral a damping term at large negative 
r Pig [15]. 

In [38], we added in a damping factor e~^'^ by hand to regulate the integrals, but this 
tends to systematically underestimate the resulting integrals, and (3 needs to be chosen 
carefully to balance accuracy and computational efficiency. However, we can regulate these 
integrals analytically, so that the numerical evaluation involves an explicitly finite integrand. 
We start by splitting the integral in Eq. (13. lip into two parts 



TO 



'^end 



h + l 



2 ■ 



(3.12) 



TO 



Here tq is an arbitrary time when all three modes are well inside the horizon, and Tend 
corresponds to a moment long after horizon exit. It is straightforward to numerically evaluate 
I21 but Ji suffers from the cut-off dependence we mentioned earlier. 

When all three modes are well inside Jhe horizon, their phase and amplitude is well- 
described by the WKB approximation 



1 



exp 



'a(k)d 



+ C.C 



(3.13) 



where a{k) = k"^ + z"/z. Deep within the horizon, k'^ ^ z"/z, the modes would not see the 
curvature term and hence it will propagate as a plane wave Vk oc exp[ikT]/^/k. Specializing 
to the case ^ = C for now 



To 



dTd{l 



y/8k^Mh 



where K = ki + k2 + k-^ and we set 



2k 



JK{t-to) 



ik 



2k 



(3.14) 



(3.15) 



-'^The 2-point function is computed with WKB methods in [48 . Here we only need the mode evolution 
inside the horizon, while a full treatment requires matching across the moment of horizon crossing, using 
standard turning point techniques. 
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Figure 1: This plot compares the 3-point correlation function, computed using the (3 de- 
scribed in [3B] and the boundary regulator described in this section, in the case of the step 
model [Sec. O with model parameters (c, d, (f)s) = (0.0018, 0.022Mp, 14.84Mp)]. We plot the 
dimensionless variable G/k^ defined in Eq. (14.201) for the equilateral case. The solid and the 
dashed hues are results obtained /? = 0.01 and /? = 0.005 respectively, while the dotted line 
is the result obtained from using the boundary regulator. The value of (3 is chosen such that 
it gives optimal results at around A; = 1; if /? is too small the early time oscillation will not 
be suppressed while if j3 is too large it will over-suppress the high k values. The boundary 
regulator does not suffer from this arbitrariness, and matches the limit found when /3 — 0. 



and 9 is some function of r given by 

e{T) = ^gie,v,v')- (3.16) 
— oo, but there is always some finite value of p such that 



In general, 9 diverges as r — 
9{t)tP — > as r ^ — oo. 

Now integrate Eq. (13.14p by parts 
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y/8kik2ks K 



9{T)e 



jK(T-ro) ro 



d9 1 -i 
dr— ■j^=^= — I 



(3.17) 



The boundary term at — oo is apparently divergent, but if we use the same e rotated contour 
as before, so — r — oo(l -|- ie) the term vanishes for any finite value of e. The remaining 
integral is now more convergent at large negative r since the integrand has picked up a factor 
oc T^^ . Integrating by parts a second time we find 



1 



y/8kik2ks 



—I 



K 



—I 



d9 



K / dr 



+ 



"° , d^9 
dr 



dr^ ^/8hk2h \ K 



\ JK{t-to) 



. (3.18) 
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The resulting integral is now convergent for any F we encounter in this work, and can be 
evaluated efficiently. We could repeat this process to further speed the convergence of the 
remaining integral, but in practice this is not necessary. We illustrate the effectiveness of 
this "boundary regulator" in Figure [H 



4 Horizon scale generation of non-Gaussianities 

In this section, we explore models of inflation where the non-Gaussianities are generated 
when the modes cross the horizon. Typically these models require a violation of slow-roll 
at some fixed physical scale. All modes experience a temporary boost in their coupling 
strengths courtesy of this violation, but only modes exiting the horizon as the violation 
occurs receive a boost in their non-Gaussian signatures. Modes deep within the horizon 
are still rapidly oscillating, and the net contribution to their non-Gaussianities cancel. The 
violation of slow roll can have two origins. The ffist is a potential with a localized feature, as 
discussed in [38], and we now generalize this analysis. Secondly, if the duration of inflation 
is such that initial transients in the dynamics have not been erased before observable modes 
leave the horizon [39] we again find a significantly boosted 3-point function at these scales, 
even though the potential is smooth. 



4.1 Features in the inflationary potential 

Consider a small step in the infiaton potential [36]. In the limit that the step is a genuine 
discontinuity, the change in the potential energy across the step would be entirely converted 
into the infiaton's kinetic energy, which is then damped away. For realistic models AV/V < 
0.01 so 0^/2 ^ V across the step, and e <^ 1. Recall that we are working in the Hubble slow 
roll formalism - if we had defined e oc [V^/V)"^ this quantity can become large across the 
step. Further, rj is the rate of change of e, so 1] and t]' can become large, provided they do so 
over a small enough range in 0. Features are thus associated with a characteristic physical 
scales and thus generate scale-dependent power spectra and higher correlation functions [SS] ■ 
The 3-point correlation function ((^(ki)^(k2)C(k3)) is dominated by the erj' term [38] . 

n«*(^-d)j j ^"'drahf]' (^ul{T)u;{T)^u;{T) + sjm^ {27if6\J2^i) + c.c. . (4.19) 

In [38], we introduced the Q to describe non-Gaussianities with both shape and scale depen- 
dence: 



(C(ki)C(k2)C(k3)), (4.20) 



kik2k3 (53(ki + ks + kg) P2(2vr)7 

where P is a constant, and for convenience we set it to be roughly equal to the magnitude 
of the power spectrum, = 4 x 10~^°. In the absence of the sharp feature, (I4.20p reduces 
to the local form with Q = (3/10) /^f' X) kf in WMAP's convention. We now construct an 
analytic approximation to this function. 

To illustrate our approach, we consider two specific features, the step [36l [38] 

1 



V^(0) = ^m202 



1 + c tanh 



d 



(4.21) 
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Figure 2: The evolution of e x 77' with units k^:T = —1.2 where fc* = 1 is set to be the scale 
when (j) crosses the center of the feature for our two models . The step (c = 0.0018, d = 
0.022Mp,(f)s = 14.84Mp) is the sohd line, and bump is the dashed line (c = 0.0005,6? = 
O.OlMp, 06 = 14.84Mp). In both cases rj' is boosted by 0(1000) boost for around one Hubble 
time {6t ^ 1 in our units). 



and the bump 



1^(0) = ImV' 



1 + c sech 



d 



(4.22) 



where c, d and (pb again respectively determines the height, width and location of the feature. 
In the latter case, c must be small enough to ensure that the field point does not get trapped in 
a local minimum. We present the numerical results for erj' in Fig. [2] and the non-Gaussianity 
profile in Figure [31 

As we see from Fig. [21 77' is non-trivial over a small range of r, but we do not have an 
analytic result for Mi(r), and in [3B] we performed these integrals numerically. Now consider 
a series of hat functions, r^' = 77^ = constant, for Tm — < t < Tm + Srm, and otherwise. 
Outside the range of these hat functions, 77 = 0(e) ~ 0, so 

f dTi = 2j2v'Jr^^0. (4.23) 

Namely, the area under the hat functions or the numerical curve of 1]' should sum to approx- 
imately zero. We approximate ({ki) by its unperturbed form 

C(A:,) = ^£=(l + zA;,r)e-^'=-, (4.24) 
although the actual mode functions can deviate from this simple form by as much as a factor 
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Figure 3: A sample result of the step (left) and bump (right) models, where for both plots 
ks = 9. We numerically computed the 3-point correlation functions of both models for 
1 < k < 9 such that k = 1/1.2 corresponds to the scale of the feature at 0<j = 14.84Mp. 
For the step model, we use (c = 0.0018, d = 0.022) while for the bump model we use 
(c = 0.0005, d= 0.01). 



of two (in addition to phase shifts) for short periods of time |36] . Using these approximations. 



(C(ki)C(k2)C(k3)) 



^ — ^ /_ r 

^ ^ Tm (J'm 



2AKt 



+sym + c.c. 



(4.26) 



Y^i^kj^i^ kik2k 



kl ImEi(iii'r) 

i 

+STm 



K 



sin Kt 



K 



sin Kt — Kt cos Kt) 



where 



K = ki + k2 + k^ 



(4.26) 



(4.27) 



and ImEi (ii^r) denotes the imaginary part of the exponential integral function. For Kt ^ 
1, ImEi(^irr) ^ -tt - cos{iKT) / {Kt) + 0{{Kt)-^)] for Kt < 1, \mEi{iKT) ^ Kt + 
OHKt)'). 

For large K, KTm ^ 1, Eq. (14.261) is dominated by the last term. 



kik2k^ 4 



Tm sin KTm sin K5Tm . 



(4.28) 



15 r 
10 
5 



-\ // » /' 

I ; / ' ^ lU- 

\^ \ 4/ ' 

\ / 



'/ M 
// x\ 



/ \ 



/ \ 
; \ 
/ \ 
I \ 

/' \i // V i/ \i 
\l /( ? 1/ 

"~k /' \\ /i ! / , 

j+v — ' /' ' \ i — ^ — ^^ — 1 — ^ — ^ — — 1 k 

2 \ s / 4 / ; \ 6 J \l 8 j '\ 

\ / ' "i /' ' / 



Figure 4: A comparison of the ansatz Eq. (]4.3ip (full line) to our numerical results (dashed 
line) for the step model (c = 0.0018, d = 0.022) on the left and the bump model (c = 
0.0005, d = 0.01) on the right with k^, ^ 1/1.2. This is a reasonable match to our analytical 
from. As explained in the text, the drop-off at small and large k is not captured by the 
ansatz but can be easily incorporated by adding in a window function when comparing the 
ansatzen with data. We have added a phase factor into the ansatz to synchronize with the 
numerical results; this phase factor is physically important but is not estimated analytically. 



Numerically, the dominant scale-dependent oscillation comes from the term sin Krm, where 
Tm is the center of the hat function, and corresponds to the moment the inflaton crosses the 
feature and 4>{Tm) ~ 4>s- We denote this scale — 1/A;*. The amplitude is further modulated by 
sin KSt^- Typically St^ <^ t^, so this terms is less important unless the non-Gaussianities 
are found over a large /c-range. From fl4.23p . we see the amplitude of (14.281) is roughly t]'At 
(At is the duration of the feature in r]'). This amplitude is consistent with the order-of- 
magnitude estimate of Ref . [3H] . 

Before we write down the ansatz, let us take a step back and consider several other limits 
not well described by Eq. (14.281) . For small K, KTm <^ 1, and the leading terms (the first 
two terms in (I4.26P ) are proportional to Ylm^m^'^"T- ~ ^'^ next order terms dominate. 



kik2k3 



(4.29) 



The non-Gaussianities vanish as Kt^ — > as these modes are outside the horizon as the 
inflaton encounters the feature. 

In the squeezed triangle limit, ki = k2 = k and fcs — 0, the second term in (I4.26P 
dominates instead of the third and we have a different behavior 



kik2k3 ^ 4/c: 



cos 2kTm sin 2A;5r„ 



(4.30) 



Comparing ( ]4.28p and ( 14.30p . we see that the latter becomes more important if one of the 
momenta /ca is smaller than l/|rm|. This is less important observationally than Eq. (14.281) . as 
it requires measurements on two widely different scales. In addition, Eq. (14.281) is incomplete 
also because it does not vanish for large and this is an artifact of the sharp edge of our 



10 



hat functions. We expect fl4.28l) to decay away after a width /\K > 1/5t where 5t is the 
smoothing of the sharp edge of the hat function. 

With these caveats in mind, we can now write down the ansatz following Eq. (14.281) by 
defining 

^^'''^^ = /feat sin + phase ) , (4.31) 



where we match the phase to our numerical results. We could improve this approach by 
taking a sum of a series of hat functions, so the final form of the 3-point correlation function 
would then be a sum over oscillatory functions 

^feat X ^ /(m) 



kik2k^ 



E /2 sin + phase) , (4.32) 

m. ^k^ J 



where fci™'' is the scale associated with the center of the m-th hat function, or further polish 
this result using an Euler-Maclaurin expansion. However, our primary goal here was to 
produce an analytic expression that mirrored the qualitative form of the 3-point function 
generated by a step, and in this we have succeeded. This expression is scale-dependent and 
thus differs from "local" and "equilateral" forms. Finally, this ansatz is clearly factorizable 
and in future work we will use it - in combination with the algorithms of [31] - to probe the 
ability of specific CMB experiments to recover this type of signal from data. 
For the step potential [38] . 

7^3/2 

fe. ~ ^ . (4,33) 

where c and d are the height and width of the step respectively. This feature is approximately 
localized, with a single scale /c*. The derivation of the phase factor is not instructive and 
we obtained it by matching to our numerical results. Refs. [501 El] use a step to explain the 
"glitch" seen in the temperature Ci for I ~ 30, and find best fit parameters (c = 0.0018, d = 
0.022, 0, = 14.84Mp). Looking at Fig. 1 we find 



^fcat 



7 sin ( -I- phase ) . (4.34) 



kik2k3 VO-83 

The amplitude is within the order of magnitude that we expected from Eq. fl4.33p and = 
— 1/k^ ^ —1.2 is indeed around the center of the feature in Fig. [2l These approximations 
are intended to capture the qualitative form of the 3-point function and do not have a high 
degree of numerical fidelity, but could certainly be further developed. Moreover, the detailed 
numerical match to the exact result could be improved substantially by adding heuristic 
parameters to the approximation, and then varying these to optimize the approximation. 



4.2 Non-attractor initial conditions 

Non-Gaussianities can also be enhanced if the infiaton trajectory is initially displaced from 
its attractor, slow- roll solution - for instance by giving it large kinetic energy ( "fast-roll" ) or 
perturbing the field point in a steep direction, orthogonal to the infiaton trajectory. Hubble 
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friction will erase this transient, so this situation only arises when the overall duration of 
inflation is close to the minimal value. Again, modes that are crossing the horizon during 
this period have their non-Gaussianities boosted due to the enhanced coupling. In this case 
the relevant scales are the longest modes that contribute to the CMB. 

This mechanism was used by Contaldi et. al. [19] to explain the possible low-/ suppression 
in CMB (see also [52l[53])| The inflaton is started with a large velocity and the potential 
is negligible, so 



+ 3H^ ^ , 




One finds a{t) ~ t^^^, and (p ~ —VQ/ (3(t + to)), where to is determined by the initial velocity. 
During this fast-roll period, 

ej ^3 , r]f ^0 . (4.36) 

In a qualitative estimate, Ref. assumes that the slow-roll inflation period begins im- 
mediately after this kinetic-energy-dominated fast-roll (non-inflationary) period, and during 
slow-roll 



0.01 . (4.37) 



Let us estimate the non-Gaussianity. For example, consider ( ]A.69I) 

■^Tbcgin \ i /V 2 / j 

(4.38) 

The fast-roll period is not inflationary, so modes are entering the horizon during this phase. 
As we will explain more quantitatively in Appendix [B], we are interested in the modes that 
are slightly within the horizon as inflation begins. Modes that are near the horizon have 
smaller non-Gaussianities. Although the coupling e/ ~ 3 is greatly enhanced comparing 
to eg ~ 0.01 in the slow-roll phase, the mode function Uk oc l/i/e/ is greatly suppressed. 
Looking at (14.381) . we find the following factors of e in the amplitude of the 3-point correlation 

^ - /fast - e}e;y\f\l . (4.39) 

The first factor is the large coupling; the second is due to the asymptotic value ^^(Tend), 
which is determined by the slow- roll inflation period; the third is from ma:(t) during the fast- 
roll period - note this is where the suppression comes from; the fourth is from the prefactor, 
1/P^ oc e^, in the definition of Q/k^. This gives 

/fast ~ 0.5 , (4.40) 

^Ref. |54| extended this approach to other non-inflationary initial conditions, and the same conclusion 
applies to this model. Conversely, [55j looks at a perturbation in a two-field model, which is beyond the 
scope of the current analysis. 
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which is confirmed by a more detailed estimate in Appendix [Bl 

We next look at the transition period from the end of the fast-roll to the slow-roll, where 
e drops from ^ 3 to ~ 0.01. This period is absent in the above analytical model of 
Ref. [H], but it also has important contribution to the non-Gaussianities. Once the kinetic 
and potential energies are roughly balanced, H is approximately constant, and (j) drops as 
g-3i?i g^j^i^ g drops as e~^^*. The r] grows from to 0{1) at the beginning of this period, 
then drops to 0(0.01) in the end. The most important contribution to the 3-point is the ef] 
term, similarly to the sharp feature case, and possess the amplitude 

Uan = OiAr]) = 0(1) . (4.41) 

The shape and running of this non-Gaussianity is also similar to those of the sharp feature, 
hence the ansatz (14.311) applies. 

The sum of these two contributions implies that the 3-point is (9(100) larger than the 
/ ~ of standard slow-roll inflation. Our numerical simulation (not shown here) confirms 
this expectation. This is only comparable to the f^j^f^^ = 0{1) 3-point expected from simple 
non- linear gravitational effects of [1], and an order of magnitude smaller than the 3-point 
we found for a bump. Since the kinetic energy is simply redshifted away by the expansion 
of the universe, we cannot easily enhance this signal by tuning the potential. Further, this 
signal peaks at scales where cosmic variance will be largest, so we are not optimistic this 
signal will be observable, even under ideal conditions. 

5 Sub-horizon generation of non-Gaussianities: Resonance model 

We now turn our attention to the generation of significant non-Gaussianities while the modes 
are still well within the horizon. Again, we inspect the integral representation of the 3- 
point correlation function, Eq. (I3.1ip . Modes that are well within the horizon oscillate 
rapidly. Since the interaction couplings e, t]' are roughly constant for plain vanilla slow-roll 
model, these oscillations cancel. However, if the interaction couplings oscillate, they can 
interfere constructively with the rest of the integrand, yielding an enhanced 3-point signal. 
These contributions are generated while the modes are deep inside the horizon, and are thus 
physically distinct from the situation explored previously. 

When all modes are inside the horizon, each of the 3-point integrals ( 13. lip consist of an 
oscillatory piece ~ e*^'^^"'"'^^"'"'^^'*'^ and a prefactor g{t,r]') which is a function of the slow roll 
parameters. If the potential has a small oscillatory component, g becomes 

^(e, T]') ~ a{l + 5sin(u;r)) , (5.42) 

where u > H, and 5 <C 1. We assume that a, uj and 5 change slowly over a single Hubble 
time, and treat them as constants in what follows. As the physical wavelength a{t)/K 
increases, the mode will briefly resonate with the coupling term g when its frequency is 
roughly uj. During this resonance phase, the ^i^2^3 term in the integrand of (13.111) can 
interfere constructively sin (cur) and generate a large 3-point function. Conversely, destructive 
interference will generate no extra contribution to the 3-point function. See Fig. O for an 
illustration of the resonance effect. 
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z 



Figure 5: This figure illustrates resonance between the total mode K = ki + k2 + momen- 
tum and an oscillatory rj'. The solid line shows C(^)) the dashed line describes 7]', and the 
dotted line is the integral Eq. (13. lip from — oo to time r with ki = k^ = k^ = k. Resonance 
occurs when the frequency of rj' is roughly 3A;. This figure was generated after numerically 
evaluating the mode functions for the parameters of Eq. (I5.46p . with all the lines rescaled to 
arbitrary units to emphasise the effect. The universe grows by roughly one e-fold over the 
range of this plot, and the relevant modes are well inside the horizon. 

This resonance requires 



in order to ensure that the relevant modes are sub-Planckian during the resonance epoch. 
In practice, this is not a strong constraint, but the resonance introduces a new length-scale 
into our analysis of cosmological perturbationsif] In this case, we expect non-Gaussianities to 
be present at all scales, since the ripple is laid down across the entire potential, in contrast 
to the "features" considered above. Moreover, the modulation of the potential must be 
small enough to ensure that the infiaton does not get trapped - and if we make it small 
enough, we can also ensure that the 2-point function is not significantly modified even if the 
non-Gaussianities are large, as discussed in Appendix O 

To study this mechanism explicitly, consider a standard slow-roll model with a very small 
oscillatory component 



Since is slowly rolling, the resonance effect occurs when the physical frequency reaches 

uj ^ — — = — — . (5.4J 

27rA V67rA ^ 

•^In models with a non-trivial UV cutoff below the Planck scale, this scale would replace Mp in (|5.43p . 



H <uj < M, 



■p ' 



(5.43) 




(5.44) 
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Therefore to satisfy the resonance condition fl5.43l) . we need A0 -C 1 and m -C A, while 
we take H = rrKp/y/G, in accordance with slow roll. In Appendix Owe show that we need 
C0/A <^ 1, so the perturbation to rj is small. 

Before we proceed, we emphasize that f l5.44p is simply a toy model, with no direct physical 
motivation. However, in brane inflation, it has been argued that the duality cascade within 
a throat will leave tiny sharp features in the inflaton potential or the background warp factor 
[56] . Such features typically come in a series and can have cosmologically observable effects 
as the inflaton branes roll across them |39j. If these features are large and well-separated, 
they would yield of sequence of isolated sharp features. However, in the opposite limit, one 
would have a situation closer to the modulated potential above. 

We now estimate the resulting 3-point function. For definiteness, we use the following 
parameters in our numerical examples 

m = 3 X lO^^Mp , c = 5 X 10"^ , A = 0.0007Mp , ^ ^ 15Mp . (5.46) 

The non-Gaussianity oscillates with K because g{e, rj') has a continually changing phase. The 
time it takes the inflaton to cross one "ripple" is At = ^/QirA/m, during which AN^ = tt^A. 
So the oscillation period of this non-Gaussianity in i^-space is 

AK = KAN^ = nK(l)A . (5.47) 

Using the parameters in fl5.46p . this gives AK/K = 0.033. This agrees with the numerical 
results shown in Fig. [61 We emphasize that the property AK oc K is model-independent, 
since, ignoring the slow-variation of background evolution, the analysis for At and AN^ is 
scale-independent. Another universal property is AK/K < 1, via fl5.43p . 

We now estimate the amplitude of these non-Gaussianities. The frequency of each mode 
is continuously decreasing due to the expansion of the universe. Once it differs from the 
resonant frequency by Acj, the integration in the 3pt starts to cancel out if it is performed 
over Ati ~ tt/Auj. Meanwhile, it takes At2 ~ Auo/{ujH) to stretch a mode sufficiently in 
order to induce a frequency change oj to uj — Auj. Equating Ati and At2 gives the time 
period over which resonance occurs for this particular mode. 



At~^/-^. (5.48) 



In f l5.44p . this corresponds to the number of oscillation cycles 

At 1 



T VA0 



(5.49) 



we need to integrate over, where T = ^/QiiA/m is the oscillation period of the background. 
For fl5.46p we need about 10 cycles, which is confirmed numerically in Fig. [51 

The dominant source of 3-point terms - as for the "feature" models - is the erj' interaction 
term Eq. fl4.19p . Because the mode is well within the horizon \Kt\ 3> 1, in 04.250 the last 
term dominates. Denoting the oscillatory behavior of the slow-roll parameter as 

V 3 ?7osci = Va sin(t^t) , (5.50) 
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where the subscript "A" denotes the oscillation amplitude. Integrating over one period 

dr T sin(cjt) e'^^ (5.51) 



gives the amplitude vrr^/iT, where = —1/aH is evaluated around the resonant point. 
Combining this with Eq. fl5.49p , from (14.251) we get the estimate of the amplitude 



To evaluate (I5.52p . we need fj. For (15.441) the amplitude of the dominant oscillating term 
in fj is (see Appendix [C]) 

/-cmd) 

Va ~ . (5.53) 

Note one might naively expect that rj ~ V" ~ c/A^, i.e. that it scales as A^^. However this 
over predicts the amplitude by a factor of A~^ : although r] oc H, since H = —0^/2 via the 
scalar field equation of motion (12.31) and hence H = = ^HcfP' + V'cj) oc V . It follows 
that f] (X H (X V". The crucial point is that one can always use Eq. (12.31) to eliminate one 
derivative of V in the derivation of H; in other words one cannot simply neglect the (p term 
in the equation of motion. Numerical computations of rj show that this estimate Eq. (I5.53P 
is accurate in both amplitude and scaling. 

Combining (I5.52p and (I5.53p . we get the amplitude 

9 cM^ 

~ 4^2.500.5 ' (5.54) 

where we have restored Mp into the final answer. For (I5.46p . /res ~ 22, which differs from 
matches to our exact numerical results shown in Fig. [2]by ~ 30%. The numerical results 
also exhibit the Q/k^ oc c/ h?'^ scaling, as can been seen in Fig. [3 
The ansatz for the non-Gaussianity is therefore 

Q 

, , , = /res sin(C In i^r + phase) , (5.55) 

1^1 1^2 1^3 ansatz 

where 

C = ^ 2/(0A) . (5.56) 

The In K comes from the fact [in Eq. (I5.47P ] that the oscillation period AK is proportional 
to K and the background dependence of AK/K is approximately scale- invariant. The 
value of C is determined by the oscillation frequency in the K-space. We also point out that 
the power spectrum has the same oscillatory frequency in i^-space, due to the same reason 
that we stated before (I5.47p . See Fig. [H]in Appendix O for details. We compare this ansatz 
with the numerical result in Fig. [61 This ansatz is factorizable at least by Taylor-expanding 
In in a range smaller than K. Note unlike the ansatz for the feature (I4.3ip where the 
non-Gaussianities are peaked around some fixed scale k^, the amplitude does not decay and 
is present at all scales. In addition, K > AK so it oscillates faster. As we noted above, the 
specific potential we consider is a toy model, but this analysis could be extended to more 
complicated modulations, or models where the modulation applied only over a finite range 
of field values. 
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Figure 6: The numerical results for the 3-point correlation function of the resonance model. 
The plot on the left hand side is a 2-D slice of the 3-D function with fixed = 1.1 while 
k2 and fcs runs from 1.0 to 1.1. The plot on the right hand side is a comparison of the 
numerical result (full line) and our analytical ansatz (dashed line) for the equilateral case 
with 1.0 < < 1.5 - we have increased the range to more show the 1/K dependence of the 
frequency. The amplitude fl5.54p over-predicts by about 30% - here we are have used instead 
numerically computed value of 15 to fit the plot better. On the other hand, the frequency 
f l5.56p is accurate - we have used C = 2.05/ {(pA) and added a phase factor to synchronize 
the plots. 

6 Discussion and Conclusion 

We studied two distinctly mechanisms for generating large non-Gaussianities within single 
field inflation. We derive the approximate 3-point correlation functions using semi-analytic 
methods, which are in the computationally useful factorizable form. In the first mechanism, 
non-Gaussianity is generated at horizon crossing by either a feature in the potential, or an 
initial transient in the inflationary dynamics. In the second case, oscillating slow roll param- 
eters induce a resonance which leads to the generation of a non-Gaussian signal well before 
horizon crossing. In both cases, the 3-point function depends strongly on the individual 
wavelengths of the modes in the "triangle". With a single, sharp feature or non-standard 
initial conditions, the 3-point is only enhanced in modes which are crossing the horizon 
as the infiaton traverses the "feature". In a resonance model such as f l5.44p . the physical 
non-Gaussianity is present at all scales. 

In the resonance case, we showed that the 3-point function is periodic with a period 
AK proportional to, and smaller than, K = ki + k2 + k^. This non-Gaussianity will peak 
starting from a fixed scale when projected onto the CMB sky. Assume for simplicity that 
K i where i is the CMB multipole, and denote Ai as the oscillation period. At larger 
scales where the oscillation spanning Ai < 1, this non-Gaussianity presumably cannot be 
resolved experimentally. For the numerical example considered here, it would become visible 
at £ ~ C(IOO) where Ai starts to exceed C(l), inducing an effective scale-dependence in 
this signal. 
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Figure 7: The numerical results for the set of parameters (c = 5 x 10 ^, A = 
0.0007Mp,shortdashed), (c = 10 x 10"^ A = 0.0007Mp, full) and (c = 10 x 10"^ A = 
0.0014Mp, longdashed). The amplitude of the 3-point correlation function is proportional 
to c/A^-^ while the frequency is independent of c and proportional to 1/A. 

With a feature in the potential, the resulting transient violation of slow-roll generically 
leads to an oscillatory and scale-dependent 3-point function. We wrote down a heuristic 
and factorizable scale-dependent expression for this signal and showed that it captured the 
qualitative properties seen in the exact numerical evaluations of the corresponding integrals. 
The 3-point correlation decays as we move away from K k^, where k^, corresponds to 
the scale leaving the horizon as the inflaton traverses the feature. We also show that while 
non-standard initial conditions such as the fast-roll model of [19] can generate large non- 
Gaussianities relative to our usual expectations for single field inflation, the amplification is 
not expected to lift them above the "noise" of non-linearities produced by post-inflationary 
gravitational evolution alone. 

In addition, we have presented detailed description of a general numerical method for 
computing the 3-point correlation functions of primordial perturbations from canonical single 
scalar field inflationary models with arbitrary potentials. We show that while the integrals 
themselves are formally convergent, they need to be regularized as the integrands are oscilla- 
tory, and show how this can be accomplished analytically, rendering the numerical integrals 
rapidly convergent. 

There is much further work to be done. Our immediate goal [57| is to use our heuristic 
ansatzen to construct an optimal estimator with which to search for scale-dependent non- 
Gaussianities in the cosmic microwave background, and estimate the likely bounds that 
future missions can put on this signal. On the theoretical front, we plan to investigate the 
details of non-Gaussianities generated by multi-field models [5[ [7] within the horizon, and 
with non-standard kinetic terms. Moreover, while the analytic approximations we present 
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here capture the quahtative form of the 3-point function, they are not intended to provide 
a precise quantitative match to the numerically computed values, but these approximations 
can certainly be improved. 
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A 3-point correlation functions 

In minimally coupled single field inflation, the cubic interaction Hamiltonian for the scalar 
perturbation C is [21 HEl fT5] 



Here d is the inverse Laplacian. This cubic Hamiltonian is exact for arbitrary e and r/. 
The 3-point correlation function at some time r after horizon exit is 



(C(r,kOC(r,k2)C(r,k3)) = -z / dr' a ([C(r, ki)C(r, k2)C(r, kg), i/,„i(r')]) , (A.59) 




(A.57) 



where 



(A.58) 




together with a term coming from the field redefinition 



(C(kOC(k2)C(k3)) 



= (Cn(kl)Cn(k2)Cn(k3)) 

+ r^(C(ki)Cn(k2)Cn(k3)) + sym + 0{r)\P^,f) , 



(A.60) 



where Cn(k) denotes the Fourier transform of 
We evaluate it using the decomposition 



C(r,k) 



u{t, k)a(k) + u*{t, 



Vk = ZUk , 




(A.61) 
(A.62) 
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where v{t, k) is the solution of the hnear equation of motion of the quadratic action, 

v'k + k^k - —Vk = . (A. 63) 



Our choice of vacuum imphes that the initial condition for the mode function is given by 
the Bunch-Davies vacuum 

Vkiro) = 

<(ro) = -z^ (A.64) 

where we have neglected an irrelevant phase. 

We get seven contributions. The terms proportional to arise from the a^e^CC^ term 

/Tend ( ^„ \ / A/* A;* \ . 

dr \\\u,{T,ni) \ (^«I^^ + two permj (27r)=^53( J] k,) + c.c. , (A.65) 

the a(?C,{dCf term 

-2i y^^T a^e^ {^Ui{T^nd)u*{T)^ (ki-k2 + twoperm)(27r)3<53(J]]ki) + c.c. , (A.66) 
and the —laeC^idQipx) term 



2i dr ah' (|n«.(re..)j (^^^^^^ + ^^6 perm) {27:f5\'£\.,) + 



The term proportional to ef] is 



c.c. . 
(A.67) 



^ ( n^^(^^"<^) I / c?™^e7;' {uI{t)uI{t)^uI{t) + two perm ) {2'Kf5^C^\^i) + c.c. 



dr 

(A.68) 

The terms proportional to include that from e/2adC,dxd''x term 



(A.69) 



and that from ^(9^C)(^x)^ term 



J^dT ah' (|n«.(re..)j [<^^kl^^ + two perm) (27r)='53(^ k,) + c.c. . 

(A. 70) 
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The field redefinition (1A.60P contributes 



V I |2| ,2 
2 l"2| \U3\ 



(27r)^5^(^ kj) + two perm 



(A.71) 



In these equations, the "two perm" stands for two other terms that are symmetric under 
permutations of the indices 1, 2 and 3. 

In Ref. pi HHl |T5] where sharp features and non-attractor initial condition are absent, 
terms flA.65p . flA.66p . flA.67p and flA.7ip give the leading contributions to non-Gaussianities. 
In Ref. [38j where features are present on an otherwise fiat potential, term (1A.68P is the 



leading term. For the non-attractor initial condition that we consider in Sec. 14. 2[ the terms 
f lA.69P and (]A.70I) can also give important contributions in addition to other terms, because 
e and t] can be large initially. 



B Non-Gaussianities from initial conditions 

In this appendix, we give some details on the estimate of non-Gaussianity from the fast-roll 
period discussed in Sec. 14.21 In terms of the conformal time r, the scale factor of the fast-roll 
period becomes 

a = ^^/l + 2hsr , Ti<T<0. (B.72) 

Here the fast-roll starts at r = Tj and end at r = 0. hg is the conformal Hubble parameter 
h = a' /a evaluated at the beginning the slow-roll inflation period (i.e. the end of the fast-roll 
period) r = 0. Hs is the corresponding Hubble parameter in terms of t, H = h/a. This 
fast-roll period is immediately connected to the inflation period in which the scale factor 
grows as 

^^/^^ 0<r<l//i, . (B.73) 



1 - hsT 

We first look at the power spectrum discussed in [19]. Using (lB.72p . one gets 

e ^ 3 , ri^O . (B.75) 

From the definition z = a-\/2e, we get 



Z (l + 2hsT) 



(B.76) 



Picking a special initial condition for Vk, the equation of motion v'j^ + k'^Vk — = can 
be solved as 
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(2) 

where Hq is the second Hankel function with index 0. Requiring this solution and its first 
derivative to continuously match to the slow-roll results 



Vk 



(r) = Cie-'''^+'''/^' ( 1 + 



-kr + k/h 

at r = determines the coefficients Ci and C2 



kr — k/h. 



[B.78) 



TC 



TT 



32K 



-ik/hs 



Jk/hs 



H, 



(2) 



(2) 



k 

k 
2h. 



+ z HI 



(1) 



k 
2hs 



[B.79) 



The power spectrum is 



JPk_ 
(27r)2e 



\Ci — C2 



(B.80) 



At large scale for k <^2hs 



H^k^ 



{2Txf thl 



In 



2K 



(B.81) 



Hence we see that ^ as ^ resulting a suppression for the large scales. For k > hg, 
it relaxes to the usual attractor value Pk = H^/{87i'^e). 

Now we estimate the order of magnitude of the contribution to 3-point function by this 
fast-roll period. As an example, we look at the terms ( 1A.69|) and ( ]A.70p . The terms 
are expected to be of a similar size. We look at the modes k > hg so the asymptotic value 
of Uk at the end of inflation Tend is approximately that of the slow-roll inflation, 



UkiTend) 



iHg 



-ikr 



fB.82) 



as we see from the discussions in the last paragraph. 

The net contribution of flA.69p and flA.TOp to the 3pt is 



64 k^/2h^J' VeJ 



3/2 



7^5^(kl + k2 + k3 



(B.83) 



where 



7^ = Re 



2/3 



dx X H^^'^\x)H'^'\x)H\ 



(2)*/ 



r(2)* 



.X 



(B.84) 



We have integrated from the start of the fast-roll period Tj to the beginning of the slow-roll 
period r = 0. We have written the lower limit of the integration in terms of the Hubble 
parameter at the beginning of the fast-roll, Hj. 
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Comparing to the WMAP's ansatz in the hmit ki = k2 = 



{^-ny^JfastPl^^X^, + k2 + ka) , (B.85) 



we get a scale dependent amphtude for this fast-roll period 

; \ 3/2 

\ a/2 3/2 



ffast ~ -0.34 (^-J ey^e7^7^ . (B.86) 

The integrand in (1B.84P starts oscillating for large x. So for very large k, k ^ 2hs {Hg/ Hf)"^^^ , 
this integration approaches zero. Therefore, we look at the modes k < 2hs {Hg/ HjY^^ . 
For example, for k/hg ^ 10, e/ ^ 3, 0.01, Hg/Hj 0.1, we have IZ = 0{1) and 

ffast ~ C(0.2). This is consistent with the rough argument given in (I4.39p . 

C Slow-roll parameters in the resonance model 

In this appendix, we work out some details on the behavior of the slow-roll parameters in 
the resonant model. Using (15.441) as the explicit example, we decompose (and similarly for 
other variables) 

(Pit) = Mt) + 0osci(t) , (C.87) 

where the (po is the unperturbed background evolution, and 0osci is the oscillation caused by 
the small high frequency ripples imposed on the potential. Since the oscillation frequency 
uj ^ H, the equation of motion for the oscillatory part is 

0osci + (V^')osci ~ , (C.88) 

where the term 3H(f)osci is much smaller than 0osci and neglected. The dominant contribution 
to (l^')osci comes from 

Ksci = ^mVcsin^ . (C.89) 

Using the background evolution 

Mt) = -^^t + (t^i ' (C-90) 

we can solve (10.880 and get 

0osci^^c0gAcos(^) . (C.91) 

With this result, one can check that flC.891) is indeed the leading oscillatory term in (15.441) . 
We next calculate the oscillatory part of H and H. Differentiating 

m'^ = V + , (C.92) 
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Figure 8: The power spectrum Pk for the resonance model. There is a small (9(10~^) 
oscillation. The frequency (in /c-space) is identical to that of the 3-point correlation function 
(in i^-space) in Eq. (15.551) . This is also different from the sharp feature case where the 
frequency (in fc-space) of the power spectrum is twice of the frequency (in i^-space) of the 
3-point correlation function [38]. 



we get 



So, we have 



From (lC.94p . we get 



H = -^02 _ (C.93) 



^ = (C.94) 
Ho 2 A 



i^osci 3c0oA fM't)^ /nnr,\ 

Note that if we try to get -ffosci from (1C.92P , the leading terms will cancel and the subleading 
terms are not reliable in our approximation. 

Now we calculate the oscillatory part of the slow-roll parameters. It is easy to see from 
the above results that the leading correction term for e = —H/H^ comes from that for H, 

so 



^ ^osci 
*^osci ~ ^0 







-3csin(^) . (C.96) 
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The leading correction term for rj = e/ {tH) comes from that for e, 



no— 

eo 



A ' A 

Taking the time derivative of the above equation flC.97|) we get 



3c0o .0o(i)x (nn^\ 
cos( — - — ) . (C.97) 



« V6!^=^si„(Ml), (C.98) 



which we can also derive directly by plugging in the potential Eq. (15.441) into the definiton 
for ?7, Eq. (12.61) . and keeping the leading oscillatory term. 

The most important term in the linear equation of motion flA.63p is 

- = 2a^H^ f 1 - £ + 3 i + i 2 il^ _ (C.99) 
z \ 2 a' A ' AH ) ^ ' 

The most significant oscillating component comes from the last term 

This term causes a dramatic oscillation of the freeze-out scale of the perturbation, which we 
define to be a^^/z/z". For example, in terms of the numerical number (15.461) . the amplitude of 
(IC.lOOp is 1.53, of the same order of the non-oscillating component 1. The horizon size that 
we used in the main text refers to the averaged non-oscillating component. This oscillation 
affects the power spectrum, which we present numerically in Fig. [HI The oscillation in the 
power spectrum is tiny, but is potentially observable if the fidelity of our CMB data is very 
high. The structure of these oscillations is reminiscent of those of the "trans-Planckian" 
models of inflation, but with the distinctive scale-dependent oscillation frequency which has 
the same properties as we discussed at the end of Sec. Ofor the non-Gaussianity. For relevant 
details on constraining such oscillatory power spectrum on the CMB sky, see Ref. 
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